Photon helicity driven electric currents in graphene 
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We report on the observation of photon helicity driven currents in graphene. The directed net 
electric current is generated in single layer graphene by circularly polarized terahertz laser radiation 
at normal as well as at oblique incidence and changes its sign upon reversing the radiation helicity. 
The phenomenological and microscopic theories of the observed photocurrents are developed. We 
demonstrate that under oblique incidence the current is caused by the circular photon drag effect 
in the interior of graphene sheet. By contrast, the effect at normal incidence stems from the sample 
edges, which reduce the symmetry and result in an asymmetric scattering of carriers driven by 
the radiation field. Besides a photon helicity dependent current we also observe photocurrents 
in response to linearly polarized radiation. The microscopic mechanisms governing this effect are 
discussed. 

PACS numbers: 73.50.Pz, 72.80.Vp, 81.05.ue, 78.67.Wj 
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I. INTRODUCTION 

Graphene, a one-atom-thick layer of graphite, was ex- 
perimentally isolated only six years ago and has since 
then revealed fascinating effects in a number of experi- 
ments owing to specifics of the electron energy spectrum^. 
The chiral motion of charge carriers leads to a pecu- 
liar modification of the quantum Hall effect^ and plays 
an important role in phase-coherent phenomena such as, 
e.g., weak localization^. The fact that the band struc- 
ture resembles the dispersion relation of a massless rel- 
ativistic particle has created enormous excitement since 
rclativistic experiments in a solid state environment be- 
came feasible^. Indeed, Klein tunneling 7 -, a relativistic 
effect predicted over 80 years ago^, has been experimen- 
tally demonstrated very recently in gated graphene struc- 
tures^. Another characteristic feature of graphene is the 
presence of two valleys, each containing a Dirac cone. 
This constitutes a two-state degree of freedom - much 
like the electron spin - which was suggested to be ap- 
plied in valleytronicsi^. Most recently several theoretical 
groups suggested that the combination of intense radi- 
ation and a constant electric field applied to single or 
multilayer graphene may result in the generation of a 
vallcy-polarizc d 12 i 13 and anomalous Hall current^. 

Here we demonstrate that the illumination of mono- 
layer graphene by radiation of a terahertz (THz) laser in 
the absence of any dc field applied to the sample causes 
directed electric currents, including those solely driven by 
the radiation helicity. Photon helicity driven currents are 
well known in semiconductor low-dimensional structures, 
and the photocurrent generation has been proven to be a 
very efficient method to study non-equilibrium processes 
in semiconductors yielding information on their symme- 
try, details of the band structure and processes of electron 
momentum, spin and energy relaxation, etcJ^— . Micro- 
scopic mechanisms of this class of phenomena in quan- 
tum wells are based on spin-orbit coupling in gyrotropic 
material o 14 i 18 ' 19 or on orbital effects originated from the 



quantum interference of optical transitions 2 ^—. Our ex- 
periments evidence that the helicity driven photocurrent 
in graphene consists of two contributions. One of them 
appears at oblique incidence only and is an odd function 
of the angle of incidence. The other has its maximum at 
normal incidence and is an even function of the incidence 
angle. We show that the first effect is caused by the circu- 
lar photon drag effect 2 ^— , see also Ref. H3, which stems 
from the simultaneous transfer of the linear and angular 
momenta of photons to the free carriers in the interior of 
the graphene sheet. The second one, however, cannot be 
attributed to any photoelectric effect in an ideal honey- 
comb lattice of graphene. This two-dimensional lattice 
possesses a center of space inversion and does not allow 
for an electric current at normal incidence of the radia- 
tion. Thus, the appearance of photocurrents at normal 
incidence is a clear manifestation of the symmetry reduc- 
tion of the system. We suggest that these currents, both 
the helicity-sensitive and those generated by linearly po- 
larized radiation, also observed in experiment, are caused 
by the edges of the real finite-size samples. 

This paper is organized as follows. In Sec. lU a short 
overview of the experimental technique is given. The ex- 
perimental results are summarized in Sec. Mil In Scc. lIVl 
we present a phenomenological description and a micro- 
scopic theory of the photon drag effect in graphene. In 
Sec. [V] we develop a theory of the edge photogalvanic 
effect responsible for the observed current at normal in- 
cidence of the radiation. 



II. EXPERIMENT 

The graphene samples were prepared from natural 
graphite using the mechanical exfoliation technique^ on 
an oxidized silicon wafer. The oxide thickness of 300 nm 
allowed to locate graphene flakes in an optical micro- 
scope and to assess their thickness. We checked the re- 
liability of this method using Raman spectroscopy and 
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FIG. 1: Photosignals, U±, in a single graphene sheet as a 
function of the angle tp, which determines the radiation helic- 
ity. The data are obtained at room temperature for various 
angles of incidence, 9q. The signal is measured nearly perpen- 
dicular to the light propagation direction applying radiation 
of the cw THz laser with the photon energy 10.5 meV, power 
~ 20 mW and a diameter of the laser spot about 1 mm. The 
data for Go = ±25° are shifted by ±25 nV [sample 1, panel 
(a)] and ±50 nV [sample 2, panel (b)], respectively. The hor- 
izontal dashed lines show a;-axes corresponding to U± = for 
the shifted data. Full lines are fits to Eq. fT}. We note that 
these fits can be obtained from superposition of the photon 
drag effect at oblique incidence [see Eqs. <|12[) . (|13|) ] and the 
photogalvanic effect at normal incidence [see Eq. I|43[l]. The 
insets show the sample geometry and experimental configura- 
tion. The ellipses on top illustrate the states of polarization 
for various angles <p. 



low-temperature quantum Hall measurements on similar 
samples 2 ^. Typically samples were p-doped by adsorbed 
contaminants in the range of n < 2 x 10 12 cm~ 2 . The 
Fermi energies were Ep < 165 meV and the mobilities 
at room temperature of the order of 2.5 x 10 3 cm 2 /Vs. 
The flakes included in this study were all single layer. 
After recording the position of the flakes with respect 
to predefined markers, we contacted them with electron 
beam lithography and thermal evaporation of 60 nm Pd 
electrodes. The resistance of graphene measured between 
various contacts was in the range of 1 kO. Some flakes 
were cut into shape using oxygen based reactive ion etch- 
ing. Insets in Figs. [1] and show the shapes and the con- 
tact geometry of the four investigated samples indicated 
by the numbers 1, 2, 3 and 4. The samples were glued 
onto holders with conductive epoxy enabling the use of 
the highly doped silicon wafer as a back gate. The sam- 
ple morphology was characterized by atomic force mi- 
croscopy measurements under ambient conditions with 
the microscope in intermittent contact mode with stan- 
dard silicon tips 2 ^. 

The photocurrents were generated at room temper- 
ature applying THz radiation of an optically pumped 
continuous- wave (cw) CH3OH laser— operating at a 
wavelength A= 118 /xm (or the corresponding photon en- 
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FIG. 2: Polarization dependences of the photoresponse ob- 
tained at normal incidence of the cw THz laser radiation for 
the samples 3 and 4. The two left panels show signals as 
a function of the radiation helicity, given by cp, and the two 
right panels as a function of the electric field orientation given 
by azimuthal angle a. Full lines are fits to Eq. (JXJ) [panels (a) 
and (c)] and Eq. ([2]) [panels (b) and (d)]. These fits can be 
obtained by the photogalvanic effect at normal incidence [see 
Eq. (J43J ] . The insets show the experimental arrangement and 
the samples geometry. The symbols on top of the left and 
right panels illustrate the states of polarization for various 
angles tp and a, respectively. 



ergy 10.5 meV) with the power P « 20 mW and a di- 
ameter of the laser spot of about 1 mm. The radiation 
was modulated at chopper frequencies in the range from 
120 to 600 Hz. The sign of the signal is defined as a 
relative phase with respect to the lock-in reference signal 
frequency, which was kept the same for all measurements. 
Additionally we used a high power pulsed NH3 laser op- 
erating at A = 148 ^m and P w 30 kW— . The photocur- 
rent induced by the cw laser is measured across a 10 MO 
load resistor and recorded with lock-in technique. The 
signal magnitudes of Figs. [T][3] are given before amplifica- 
tion. The photoresponse to pulsed radiation is measured 
across a 50 resistor. 

The samples were excited both at normal and oblique 
incidence. In the case of oblique excitation, the angle #0 
between the light propagation direction and the sample 
normal was varied from —25° to ±25°. The experimental 
geometries are illustrated in the insets in Figs. Q] and [2] 
for oblique and normal incidence, respectively. In order 
to vary the radiation helicity we used a A/4 plate. The 
rotation of the plate resulted in the change of the de- 
gree of circular polarization after P c ; rc = sin 2ip, where 
tp is the angle between the initial polarization vector of 
the laser light E and the c-axis of the plate. The light 
polarization states for certain characteristic angles ip are 
sketched on top of Figs. [T] and [2j In some experiments 
we also used linearly polarized radiation. In this case 
the plane of polarization of the radiation was rotated by 
A/ 2 plates. This enabled us to vary the azimuthal an- 
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gle a from 0° to 180° covering all possible orientations 
of the electric field vector in the graphene plane. The 
relative positions of the light polarization with respect 
to the graphene sample edges and contacts are shown on 
the top of Fig.H». 

III. EXPERIMENTAL RESULTS 

Irradiating graphene samples at normal as well as at 
oblique incidence we observed photocurrents in all inves- 
tigated samples. The signal was detected at room tem- 
perature for both low power excitation level of cw laser 
and high power pulsed laser. In the latter case the width 
of the photocurrent pulses was about 100 ns. which cor- 
responds to the THz laser pulse duration. In order to 
prove that the signal stems from the graphene flakes and 
not, e.g., from the substrate, we removed the graphene 
layer from one of the samples and observed that the sig- 
nal disappeared. 

FigureQ] shows the photoresponse U± , which is pro- 
portional to the photocurrent j, of two samples 1 and 2 
as a function of the angle tp, assigning the helicity. The 
signals are measured at different angles of incidence 
6*o = 0, ±25°. While the figure shows the photoresponse 
obtained from the pairs of contacts oriented almost per- 
pendicular to the light propagation plane, a photocur- 
rent has also been observed at any other pair of contacts. 
Generally, the helicity dependence of the photocurrent 
can be well fitted by 

J = A sin 2ip + B sin 4(p + C cos 4ip + D , (1) 

where A, B, C and D are fitting parameters. Such be- 
havior, phenomenologically well described by symmetry 
arguments (see below), was found in all graphene sam- 
ples. The fits to experimental data shown by solid lines 
in Fig. Q] demonstrate a good agreement. The first term 
on the right hand-side of Eq. ([T]), which is proportional 
to sin 2tp and described by the parameter A, changes its 
sign upon reversing the photon helicity (marked by ar- 
rows). The analysis of the experimental data gives an 
evidence that a substantial contribution to the total pho- 
tocurrent changes its sign upon switching the radiation 
helicity from right- (u+) to left-handed (u-) circularly 
polarized light, i.e., for the angles ip = 45° and ip = 135°, 
respectively. This is most spectacularly seen in the data 
obtained on sample 1 , where the changing sign of the cur- 
rent is directly detected at normal incidence, 6 n = 0° 
(Fig. [TJl) . The helicity driven photocurrent gives a sub- 
stantial contribution to the photoresponse of two other 
samples as demonstrated in Figs.[2ja) and (c) for normal 
incidence of radiation. 

A deeper insight into the photon-helicity driven pho- 
tocurrents is given by their dependences on the inci- 
dence angle 9o presented in Fig. [3] for the samples 1 and 
2. We obtained that the photoresponse of the helicity 
driven electric current is substantially different for the 
geometries, where the photocurrent is measured in the 
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FIG. 3: Signals due to circularly polarized radiation Ucpc = 
[U(<p = 45°) — U(ip — 135°)]/2 measured as a function of the 
incidence angle 8o in response to the cw THz laser radiation. 
Triangles and full circles show the transversal response, U±, 
for samples 1 and 2, respectively. Open circles show the signal 
obtained in the longitudinal geometry, U\\, for sample 2. The 
inset illustrates the experimental configuration. 

direction transversal and almost longitudinal to the light 
propagation^. While the photocurrent measured in the 
longitudinal geometry is almost independent on the angle 
of incidence and reflects the signals obtained at normal 
incidence, in the transverse geometry it shows a super- 
position of two contributions being even and odd in the 
incidence angle. An important feature of the latter con- 
tribution is that its direction is determined by the light 
propagation plane, the radiation helicity and the sign of 
(9o, irrespective to the orientation of the sample edges by 
contrast to the photocurrents being even functions of the 
angle 8 . 

Now we turn to the photocurrent contributions pro- 
portional to the coefficients B and C in Eq. ((T|). This 
photocurrent, in fact, does not require the radiation he- 
licity and reflects the linear polarization component of 
the incident radiation. This has been checked in an in- 
dependent experiment. By rotating now a A/2 plate we 
varied the relative position between the plane of the ra- 
diation polarization and the axes of the samples charac- 
terized by the azimuthal angle a. The data are shown in 
Figs.[2Jb) and (d) and can be well described in agreement 
with Eq. (P) by 

J = 2B sin 2a + 2C cos 2a + D-C. (2) 

Corresponding fits are plotted by the solid lines in 
Figs.[^b) and (d). Such a behavior of the photocurrent 
clearly demonstrates the linear-polarization sensitivity of 
the effects. Note that the difference D — C actually con- 
stitutes the polarization independent effect, which can 
be related to e.g. contacts and is outside of the scope of 
the present paper. Studying the dependence of the pho- 
tocurrent on the angle of incidence for linearly polarized 
excitation we observed that, in contrast to the excitation 
with circularly polarized light, the measured photocur- 
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rents both in longitudinal and transverse geometries al- 
ways result from the superposition of the contribution 
weakly dependent on 9q, which stems from normal inci- 
dence, and the one being odd in the angle 9$ (not shown). 
Moreover, these contributions odd in 9q are intercon- 
nected: While the longitudinal one behaves as cos 2a, 
the transverse component is proportional to sin 2a. 

To summarize, we observed that both the helicity de- 
pendent and helicity independent photocurrcnts have two 
contributions: (i) an even function of the incidence an- 
gle #0i which dominates the signal at normal incidence, 
and (ii) a odd in the incidence angle contribution, which 
appears at oblique incidence only. Below in Sec. IIVI we 
demonstrate that the contributions to the photocurrent 
odd in the angle 9q, both circular and linear, stem from 
the interior of graphene and are caused by the photon 
drag effect, where the transfer of the momentum of the 
photons results in a directed motion of electrons. We 
present the phenomenological approach, which is based 
on symmetry arguments and is free from details of the 
microscopic model, and then the microscopic theory is 
put forward. The photoresponse at normal incidence, by 
contrast, can not be related with an ideal bulk material. 
In Sec. |Vl we show that these effects are due to the sam- 
ple edges and give the corresponding microscopic theory. 



IV. THEORY OF PHOTON DRAG 

A. Photocurrents as a second-order response to 
the electromagnetic fields 

In general, the photocurrents under study can be re- 
garded as a second-order nonlinear steady-state response 



if (0,0) = a^,q)E v {q,w)E;{q^) (3) 
to the electric field of a plane electromagnetic wave 



(2) 



E(r,t) = E(q, w)c- iut+icir + E*(q,L 



^iujt—iqr 



(4) 



Here E(ui,q) is the complex electric- field amplitude with 
the frequency to and the wave vector q, ji( 2 )(0,0) is the 
photocurrent density, zeros (0,0) indicate that we con- 
sider an electric current averaged in time and space. In 
a two-dimensional system the index A runs only over the 
two in-plane coordinates x and y while the indices v and 
r\ can include the normal coordinate z. The second-order 
conductivity ct^ 2 \lo, q) can be expanded in powers of q, 
hereafter we retain only the first two terms 



<i> 



XfjLwq 



(5) 



(2) 

The tensor crj^^aj, 0) requires a lack of the inversion 
center in the symmetry point group of the system, it 
can conveniently be decomposed into two tensors <r( 2,± ) 
symmetrical and antisymmetrical with respect to the in- 

(2 -M 

dex interchange v <-> rj. The symmetrical tensor <j\' 



describes the linear photogalvanic effect (PGE). The an- 

(2 — ) 

tisymmetric tensor c^' is dual to a pseudotensor 7,^ 
describing the circular PGE^ ~ 17 i 31 . The contributions 
described by the fourth-rank tensor <fr are allowed both 
in centrosymmetric and noncentrosymmetric media. In a 
simplified qualitative picture the linear-q contribution in 
Eq. ((5|) arises due to transfer of momenta from photons 
to free carriers and, therefore, is called the photon drag 
effect. 

The photon-drag current can equivalently be presented 
in terms of spatial derivatives of the electric field E(r,ui) 
as follows 



.(drag) 



E v (r,cj) 



(0) 



dx. 



E*(r,u>) 



dE v {r,u) 



dx b 



For the plane wave Q the bilinear expression in the 
square brackets is independent of r. Equation ([5]) per- 
mits one to make sure that the tensor $A^r f in Eq. (O 
describes the drag current even in case of the sample di- 
mension L being smaller than the light wavelength pro- 
vided L by far exceeds the carrier free path length I. 

Similarly to the second-harmonic generation linear in 
the light wave vector, there are two mechanisms of the 
photon drag effect. The linear- q terms can appear in 
Eq. ([5]) cither due to the gradient of electric field (qE 2 - 
mechanism) or due to a combined action of the electric 
and magnetic field of the electromagnetic wave (SB- 
mechanism or the high-frequency Hall effect). In some 
sense, these mechanisms are analogues, respectively, of 
the magnetic-dipole and electric-quadrupolc contribu- 
tions to second-harmonic generation described by the dif- 
ference j (2) (2w, 2q) - j (2) (2uj, 0). 



B. Phenomenological description of the drag 
current in graphene 

The ideal honeycomb lattice of graphene is described 
by the point group Dq^. The group contains the space 
inversion and, therefore, allows the photon-drag currents 
only because the circular and linear PGE are excluded by 
symmetry arguments^. Decomposing the tensor ^\^ n 
into symmetrical and antisymmetrical parts concerning 
the indices v and 77 the drag current can can be written 
phcnomcnologically also as 



jx ^Xfii^r/Qfi' 



(7) 



Here / is the light intensity (energy flux through unit 
surface) related with the electric field by \E(q,u>)\ 2 = 
(2n I ' cn^l , n^, is the refractive index at the frequency w, 
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c^uri is the unit antisymmetric tensor, j is the current 
density, e is the unit polarization vector of the electro- 
magnetic wave, e is the unit vector pointing in the light 
propagation direction, P C i rc is circular polarization degree 
related to the vectors e and e by P C ii C e = i(e x e*). 

The forth-rank tensor T\^ vri describes the linear pho- 
ton drag effect which is insensitive to the sign of circular 
polarization and reaches its maximum for linearly polar- 
ized light, while the third rank pseudotensor T\^ v stands 
for the circular photon drag current which changes its 
sign upon the reversal of photon helicity. 

The symmetry analysis shows that, in the Dq^ point 
group, the tensor Tx^n has four linearly independent 
components. T± = T X xxx Pxxyy? T2 — T X xxx Pxxyy = 
2T xyX y, T 3 = 2T XZXZ and T 4 = T xxzz . Then the 
in the right-hand side of Eq. ([7]) is reduced to 



\^x\ 2 + \e y \ 2 
2 

lej 2 - |e, 



T 2 q x 



T 3 q z 



-I + T 4 q x \e z \ 2 I . 



j y =Txq y - 



+ T 2 [q y 



T 3 q z ey6z I + T 4 q y \e z \ 2 I , (8b) 



where x and y are the axes in the graphene plane, and z 
is the structure normal. Equations (JSJ clearly show that 
the photon drag current can be induced only at oblique 
incidence of the radiation and vanishes under normal in- 
cidence when both the in-plane component of the photon 
wave vector q and the out-of-plane component of the po- 
larization vector e are zeros. The photocurrent contains 
both the longitudinal component induced in the incidence 
plane and determined by the parameters T± , T 2 , T 4 and 
the transversal component proportional to T 2 . The cur- 
rent described by the parameter T 3 requires z-component 
of the photon wave vector. 

For an elliptically polarized radiation, in particular, for 
circular polarization, the photon drag current in struc- 
tures of the £>6h symmetry contains additional contribu- 
tions proportional to components of the tensor: 



Jx 



■hi 



T±qy Pqi^qC Z I 



T2q z PciYcS y I ■ 



-Tiq x P cnc e z I + T2q z P c i rc exI ■ 



(9a) 



(9b) 



and sensitive to the radiation helicity. Here, T\ = T xyz 
and T 2 = T yzx are linearly independent parameters. Sim- 
ilarly to the linear photon drag effect, the circular pho- 
tocurrent §§§ can be induced at oblique incidence only. 



However, in contrast to the former, the circular photocur- 
rent always flows perpendicularly to the light incidence 
plane. 

It follows from very general considerations that, in two- 
dimensional systems, the contributions from T 3 ,T 4 and 
T 2 are small and unlikely to be observable. Indeed, in 
terms of the semiclassical Boltzmann equation the two- 
dimensional carriers are unaffected by the normal compo- 
nent of the electric field E z , provided that the graphene 
layer is flat, and the terms in Eqs. (|8"]l. ((9]) proportional 
to e z ,e* or \e z \ 2 vanish. For intcrband quantum optical 
transitions, the photocurrent governed by the component 
q z related to the currents tx q x , y has a small parameter 
ka, where k is the electron wave vector referred to the 
band extremum point and a is the well thickness in the 
semiconductor quantum wells or the carbon atom radius 
in the case of graphene. For intraband indirect optical 
absorption involving virtual transtions via other bands 
the parameter ka is multiplied by another small parame- 
ter fiw/A, where A is the energy distance to other bands. 
In graphene A is the energy difference between the tt and 
a orbitals and has an order of several electronvolts^ 
Therefore, in the following we will ignore the T 3 ,T 4 and 
T 2 terms in Eqs. ©, ©. 

Most studies on photon drag effect have been carried 
out in crystals with cubic symmetry^—, simple metals^ 
and in atomic gase o 40 ' 41 . In these systems the circular 
photon drag is forbidden because the symmetry equal- 
izes the coefficients T\ and Tq in Eq. © so that the pho- 
tocurrent j becomes proportional to the vector product 
q x e = 0. In systems of uniaxial symmetry, anisotropic 
crystals, quantum-well structures, supcrlattices etc., the 
coefficients T± and T 2 are linearly independent and allow 
the circular drag photocurrent. Even more so, this effect 
is allowed in graphene where, as mentioned above, the 
coefficient T 2 is negligible. It is worth mentioning that 
the appearance of the transverse linear photon drag ef- 
fect in the vicinity of the metal surface was discussed in 

Refs. mm 

Figure [4] schematically illustrates the relative orienta- 
tion of all remaining contributions in respect to the light 
propagation direction and its polarization state. This 
figure being in correspondence to Eq. © shows that the 
helicity driven photocurrent appears at oblique incidence 
only and flows perpendicularly to the light propagation 
direction, see Fig. [4ja) . Furthermore it changes its sign 
upon reversal of the incidence angle. This is in full agree- 
ment with the experiment, see Fig.[3l where the odd con- 
tribution to the circular photocurrent is indeed observed 
for the transverse geometry only. By contrast, the linear 
photon drag effect is generated in both longitudinal [see 
Fig- Sib)] and transversal geometries [see Fig.0Jc)] being 
also in agreement with the experimental findings. 

Also, the polarization dependences given by Eqs. ©, 
([9]) are observed experimentally. In order to simplify 
the comparison with the experiment we rewrite Eqs. ([5]) 
and ^ for the geometry relevant to the experiment with 
the light incidence plane being (xz) and the polarization 
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FIG. 4: Schematic illustration of the photon drag effects. 
Panel (a): circular photon drag effect, where the current 
is generated transverse to the light incidence plane (see 
Eqs. @). Panels (b) and (c): linear photon drag effect. Here 
both longitudinal and transverse current can be generated [see 
Eqs. ©]. 

state of the light described by the angle ip for elliptically 
polarized light (in particularly, circularly polarized light) 
and the angle a for linearly polarized light. In the first 
case one obtains 

j x ((p,e ) oc Tlb ™ 00 [4 C0S 2 6> + sin 2 O (1 - cos4</>)] 



[4 cos 2 # + (1 + cos 2 6» )(cos \<p - 1)] 



j y (<p,6o) oc Ti sin 9 cos 9 a sin 2ip + — sin 9 cos 9 sin 4(p, 

(10) 

while in the second case one has 

j x (a,9o) oc T\ sin 9q (1 — sin 2 9q cos 2 a) 

+T2 sin 9o (cos 2a — cos 2 a sin 2 #0) 

j y (a,9o) oc sin 0o cos9qT2 sin 2a. (11) 

For the small incidence angles as in our experiment one 
retains only linear in 9q terms and obtains from Eqs. (|10[) . 
(fTT]) the following simplified relations describing the pho- 
tocurrent dependences on ip: 

j x (<p,6 ) oc Ti sin 6<o + T2hl ^ 6 ° (cos4y> + 1) 



j y (ip, 9q) oc Ti sin #0 sin 2y + — sin6*o sin4</?, (12) 



and a: 

j x (a, 9q) oc Ti sin 9q + T2 sin 9q cos 2a 

j y (a,9o) oc sin^o cos6*oT2 sin2a. (13) 

The dependence of the experimentally observed pho- 
tocurrents odd in the angle incidence on the parameters 
a and ip characterizing the polarization states of the ra- 
diation is in agreement with Eq. (fT2"j) and (fl~3"|) . Indeed, 
as it is seen from Fig. [1] (also from Fig. [3} the odd in 
9q contribution to helicity driven photocurrent, oc sin 2^, 
is observed only in the transversal with respect to the 
light incidence plane geometry. By contrast, the contri- 
bution to the current odd in the incidence angle, which 
was measured in the longitudinal geometry, is described 
by cos 4ip or cos 2a (not shown) demonstrating the linear 
photon drag effect. 

It is worth mentioning that the coefficients describing 
the linear and circular drag effects have different prop- 
erties under time reversal, t — > —t. Namely, the circular 
polarization sign changes at time reversal while the linear 
polarization sign remains. Both the electric current and 
the photon wave vector are odd functions with respect 
to the replacement t — > —t, therefore the constants T\ 
and T2 are proportional to an odd number of dissipativc 
parameters and the constants Ti (i = 1, . . . , 4) to an even 
number of dissipativc parameters. 

The above analysis, which yields a good agreement 
with the experiment, was performed using the ideal sym- 
metry for the graphene layer characterized by the point 
group D 6 h- The real structures, however, are deposited 
on a substrate, which removes the equivalence of the z 
and — z directions and reduces the symmetry to the Cq v 
point group. Our analysis demonstrates that the photon 
drag effect does not change qualitatively with the reduc- 
tion of symmetry and only the values of the parameters 
Ti,...,T4, Ti, and T2 may change. However, now the 
system lacks an inversion center, and photogalvanic ef- 
fects become possible. The photogalvanic contributions 
to the total electric current are phenomenologically given 
by 

jx = X\»v 2 1 + 7A M -Pcir C e M / • (14) 

In structures of the Cq v point group, non-zero compo- 
nents of the tensors X\^v and 7a^ describing the lin- 
ear and circular photogalvanic effects, respectively, are 
Xxxz = Xyyz and j xy = — 7 VX . It means that the linear 
photogalvanic effect requires a z-component of the elec- 
tric field. Moreover, the circular photogalvanic effect also 
needs this component because, e.g., e x ~ \{e x e* z — e z e*) 
vanishes if the z component of the field is absent. These 
requirements result in vanishingly small contribution of 
the photogalvanic effects to the total electric current, 
similar to the arguments, which exclude photon drag ef- 
fect components given by the constants T3, T4 and T2. 
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To summarize, the analysis of the photocurrents odd- 
in 9 can be reduced to the photon drag effects solely, 
and can be applied to the experimental data independent 
of the influence of the substrate indicated above for the 
symmetry of the graphene flakes. 

C. Microscopic theory 

Now we turn to the microscopic theory of both lin- 
ear and circular photon drag effect. In the following we 
focus on the classical regime of interaction between the 
radiation and the electron ensemble in graphene which is 
realized as soon as frui <C £ , where £ is the characteris- 
tic electron energy, the Fermi energy in the degenerate 
electron gas or thermal energy ksT for non-degenerate 
electrons (fee is the Boltzmann constant). In this case 
only intraband transitions are involved. This approach 
is similar to the classical consideration of the drag effect 
in semiconductors^—. The theory for the drag effect 
in the quantum range of frequencies, hio ~ £ or higher, 
involving both interband and intraband optical transi- 
tions is out of scope of this paper and will be presented 
elsewhere flij ]. 

1. Qualitative picture 

First we present a qualitative microscopic picture of 
the drag effect exclusively based on Newton's second law 
of motion 

^L + l = eE(r,t) + -[vxB(r,t)], (15) 
at t c 

where e = — | e | is the electron charge, p and v are 
the electron momentum and velocity, p/r is the friction 
force due to electron scattering with r being the scat- 
tering time. Assuming for simplicity the propagation of 
linearly-polarized radiation along the graphene (grazing 
incidence) with in-plane wave vector and electric field 
and out-of-plane magnetic field: q || x, E \\ y, B \\ z. 
as shown in Fig. [SJ The generation of a dc current can 
be described in the framework of the high-frequency Hall 
effect^. It is seen from Fig.[5ta) that during the first half 
of the period T w = 2ir/u> of the electro-magnetic field os- 
cillation where the both field components E y and B z are 
positive the electron is exposed to action of the two forces 
F y = eE y and F x = ev y B z like in the Hall effect. In the 
absence of friction the phase shift between oscillations of 
the velocity v y and the field B z equals 90° and the aver- 
age value of F x vanishes. Allowance for the friction re- 
sults in a drift along the x axis. In the second half-period 
T^/2 both fields E y and B z simultaneously reverse and 
the drift velocity retains its direction. As a result, a non- 
vanishing time-averaged drag current is induced in the di- 
rection of q. Obviously, the reversal q — > q results in a 
change of relative sign of the electric- and magnetic-field 
components and, consequently, in the current reversal. 



This is the well-known longitudinal linear photon drag 
effect. In graphene, as distinct from bulk cubic crystal, 
the photon-drag mechanism due to the high-frequency 
Hall effect has a characteristic polarization dependence, 
see Eqs. © for T 3 = T 4 = 0. Indeed, if the electric field 
of the grazing-incidcncc radiation is rotated around q by 
an angle a, then its in-plane component varies oc cos a, 
the out-of-plane component of the magnetic field B z be- 
haves just as E y , and the drag current j x oc e 2 decreases 
by cos 2 a = (1 + cos 2a) /2. The geometry of incidence in 
the plane (x, z) at arbitrary oblique angle 9q is analyzed 
in the same way. For the s-polarized light, the considera- 
tion is completely the same, a decrease of B z by a factor 
of sin #o is automatically reproduced by the correspond- 
ing change of q x = |q|sin(9o- For the linear polariza- 
tion with the angle a different from integer multiples of 
90°, the both components E x and E y are non-zero which 
yields a y-component of the drag current proportional to 
sin #o sin 2a. 

For the i?i?-mechanism, we could ignore the spatial 
dependence of the electric and magnetic fields. While 
considering the gi? 2 -mechanism we can ignore the mag- 
netic field in the right-hand side of Eq. (|15|) but take into 
account the first-order spatial correction and present the 
electric field of the linearly polarized light as 

E(r,u) = 2E [cos (ut) — 2q x xsm(ut)] , (16) 

where Eq is the electric field amplitude assumed, without 
loss of generality, to be real. At the first step we find the 
linear response 

/ -iwt \ 

pit) = er^onRe — (17) 

\ 1 — IU!T J 

and the corresponding velocity v(t), where ^o|| is the 
in-plane component of the vector Eq. At the next 
step we calculate the oscillation of the electron position 
x(t) = J v x (t)dt. Finally, the function x(t) is substituted 
into the second term of Eq. P^)l and and then the time- 
average of this term is taken. As a result we obtain an 
additional contribution oc q x E$ to the continuous force 
acting on the electron. 

The circular photon drag effect 2 ^— can be described in 
a similar way. However, in this case one should take into 
account in Eq. (TT5)) not only the Lorentz force or the spa- 
tial gradient of the electric field but also the phase shift 
by 90 degrees of the two orthogonal oscillating electric- 
field components E\ _L E2 -L q. Then one can show that, 
since the phase of the oscillation p(t) is retarded in time 
by arctan(o;r) relative to E(t), see Eq. (|17p . the circu- 
larly polarized electromagnetic wave as well induces the 
drag current which is perpendicular to qy . The sign re- 
versal of radiation helicity means a change 90° — > —90° 
of the phase shift between the E\ and E2 components 
and a reversal of the circular drag current. 

The performed qualitative consideration clearly con- 
firms the statement made in Sect. IV. A that the inter- 
pretation of photon drag current in study is independent 
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FIG. 5: Schematic illustration of the high-frequency Hall ef- 
fect resulting in the photon drag effect. 



of the relation between the sample linear in-plane dimen- 
sion, L, and the light wave length A but requires but re- 
quires a small ratio l/L, where I is the mean free path 
length vt, as well small product ql. 



2. Boltzmann equation treatment 

We develop a microscopic kinetic theory of the photon 
drag effect for the simplest model of the graphene band 
structure. We consider only the conduction and valence 
band states in the K and K' valleys formed from the 
7r-orbitals of carbon atoms. In each valley, the effective 
Hamiltonian is a 2x2 matrix 



U = hv(k ■ er) 



(18) 



which describes both the conduction and valence bands 
of a graphene layer. Here v is the electron speed in 
graphene, k = (k x ,k y ) is the two-dimensional wave 
vector referred to the Dirac point, a is a matrix two- 
dimensional vector with components a x and a y being the 
Pauli matrices in the pseudospin space. 

In accordance with the definition (fj| the electric and 
magnetic fields of the incident wave are taken in the form 



E = E c 



\{qr—u)t) 



c.c. , B = Bqg 



'i(qr-ujt) 



(19) 



where Eq and Bq are the complex amplitudes which are 
orthogonal to each other and to the wave vector q. In 
the case of a linearly polarized light the vectors q, E 
and B form a right-hand triple. It is worth noting that 
the electric and magnetic fields acting on an electron in 
graphene lying on a substrate arc different from those 
of an incident wave owing to the presence of substrate. 



To be specific, let us assume that (xz) is the radiation 
incidence plane. Therefore, E Qty = Eq (1 + r s ) and 

Eq, x — Eq\.(1 + r p ), where the superscript (i) denotes 
the incident waves and r s , r p are the amplitude reflec- 
tion coefficients in the s and p polarizations, respectively. 
Similar relation holds for the z-component of the mag- 
netic field, Bq z = Bq z (1 + r s ). We expect that the 
metallic contacts attached to our small samples cause an 
inessential distortion of the plane-wave character of the 
electromagnetic field. 

Following Ref. H3 the kinetic equation for the electron 
distribution function f(k,r,t) in a given valley reads 



dt 



dr 



-\v x B] 

c 



(20) 



Here c is the light speed in vacuum, v — is the electron 
velocity in the state with the wave vector fe, 



1 de k 



k 

W\ 



(21) 



£fe = vhk is the electron dispersion, and Q{f} is the 
collision integral. 

The drag current is proportional to the light intensity 
and, therefore, appears in the second order in the elec- 
tromagnetic fields. Correspondingly, we solve Eq. ([2TJ)) 
by iterations with respect to the electric and magnetic 
fields and represent the electron distribution function as 



f(k,r,t) = f (e k ) 

+ \f 1 (k)e i ^ r - u ^+ c.c] +f 2 (k) + 



(22) 



where fo(sk) is the equilibrium distribution function, 
fi(k) describes the linear response to the fields, f-2,{k) 
is a homogeneous time- independent correction which ap- 
pears in the second order in E, B, and the omitted terms 
(...) describe other contributions including those oscillat- 
ing with a double frequency and terms of the higher order 
in E,B. The direct current caused by the photon drag 
effect is given by 



3 = 4e^u fc / 2 (fc) 



(23) 



where the factor 4 accounts for the spin and valley de- 
generacy. 

The first-order correction fi(k) can be found by solv- 
ing the linearized Eq. (j2T))) with the result 



h(k) = - 



(ev) - ir 2 



\UJT\ 

(qv)(ev) - v 2 (qe)/2 

1 — 1UJT2 



(qe)v 2 



2lc 



(24) 



Here, /q = dfo/dsk, t± is the momentum relaxation time 
(more precisely, relaxation time of the first angular har- 
monic of the distribution function) and t 2 is the relax- 
ation time of the second angular harmonic. While deriv- 
ing Eq. (|24[) . we took into account that T2\q\v, \q\v/ui <C 
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FIG. 6: Theoretical frequency dependences of (a) the circular 
photon-drag current, and (b) the linear photon-drag cur- 
rents, both longitudinal (j'y) and transversal (j±), induced in 
bulk graphene sheet at oblique incidence and calculated after 
Eq. (|26[) . (|27[) . The drag current is given in arbitrary units 
and the frequency is presented in the dimensionless form as 
u)T\, where n is the momentum relaxation time. 



1 , in which case the second and third terms in the square 
brackets are small corrections to the first one. More- 
over, the energy relaxation time r e (relaxation time for 
the zeroth angular harmonics of the nonequilibrium dis- 
tribution function) was assumed to be much longer than 
all other time scales in the system, namely, ti,T2,1/lj. 
The latter assumption allows us to neglect corrections 
to the distribution function related to energy relaxation 
processes and allow for this relaxation time to be infi- 
nite^. In this connection we note that the first term in 
Eq. contains the first-order angular harmonics of the 
wave vector, the second term contains the second-order 
harmonics while the third term is angle independent. 

The second-order correction /2(h) averaged in time 
and space is found from the equation 



2e„ 
—Re 

h 



En 



-\vx B ] 

c 



(25) 



where the asterisk denotes the complex conjugation. 
There are two contributions to /2(fe) related to the above- 
mentioned EB- and g£? 2 -mechanisms: the first arises 
from a product of the Lorentz force and the main term 
in Eq. (f2"4")l . and the second results from a product of the 
electric-field force and the ^-dependent terms in Eq. (f2"4"f . 
Substituting the solution of Eq. (|2"5j) into Eq. ([25)1 we 



arrive at the following expressions 
47reV ^ ri/£ 



d,T\ 
de k 



To 



3„,4 



+ uj 2 t 2 

1 - W 2 TlT 2 

1 + uj 2 t 2 
' ri.fi 

^ 1 + UJ 2 T 2 



( dn 


--)} 


\de k 


£ kj . 


/dn 




\de k 


£kj 



(26a) 



(26b) 



for the constants T\ and T 2 describing the linear photon- 
drag current in the phenomcnological equations ([5]). It is 
worth noting that in the low frequency limit, wt\,u>T2 "C 
1, the coefficients Ti, T2 tend to infinity as uj^ 1 . However 
the drag current has a finite limit at uj — > since it is 
proportional to gTi j2 and q ex uj. It should also be pointed 
out that, according to the definition (fT9|) . the static fields 
E and B differ by factors of 2 from the amplitudes E 
and B of the corresponding fields taken at uj — > and, 
therefore, the constant relating the current at uj — » with 
the squared amplitude E 2 and the constant relating the 
static Hall current with the product of static fields E and 
B differ by a factor of 4. 

In the same way we derive the constant Ti describing 
in Eqs. (|9]) the circular photon-drag effect 



Ti = 



27reV v r 2 (l + T 2 /r 1 )f> 
c A.[ 1 + ( wri ) 2 ][l + (^r 2 ) 2 ] 



G?Tl 

cfefc 



n 



(27) 

It is finite at uj — > which means the vanishing of the 
circular photocurrent in the static-field limit as expected 
because, for static fields, ellipticity is forbidden. 

The time inversion symmetry imposes the certain re- 
strictions on the coefficients in Eqs. © and micro- 
scopic expressions for Tj (j = 1...4) and Tj (j = 1,2) 
must be proportional, respectively, to even and odd num- 
ber of dissipative constants, in our case the relaxation 
times T\, T2 or the inverse times t^ 1 , tv^ 1 . One can see 
that Eqs. (|2"6"|) and (|2~T|) for the linear and circular drag 
effects satisfy these general rules. 

Thus, we have derived microscopic expressions for the 
phcnomenological constants in Eqs. ([8a|). (|8b|) . Eqs. (|9aj) 
and (|9bp which, as we addressed above, describe well all 
experimental findings. The developed microscopic theory 
yields all allowed contributions to the photon drag effect 
in the model where only 7r-orbitals of carbon atoms are 
taken into account. 

Figure [B] shows calculated frequency dependences of 
the photocurrent. In the calculation we assume that the 
electron scattering in graphene is short-range, in which 
case one has T\ = 2t 2 cx eZ . It is seen that at small fre- 
quencies, oj — > 0, the circular photon drag effect vanishes 
and the linear photocurrent reaches its maximum values. 
The absolute value of j c has a maximum at ujt\ ks 0.8 and 
decreases as uj~ a at high frequencies, ujti,ujt2 1. The 
linear photon-drag current exhibits a decrease at high- 
frequency oc UJ~ 3 . 
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FIG. 7: Schematic illustration of the edge current generation 
in a sample of finite size. 



Finally, we note that together with electrons in the 
conduction band, holes in the valence band also con- 
tribute to the photocurrent. The hole contribution to the 
current is obtained from Eqs. (|26|) . ([27|) by the replace- 
ment e — > — e and the equilibrium electron distribution 



(hole) 



function /o by the hole distribution /g 



V. EDGE PHOTOCURRENTS 

Here we develop a theory of the photocurrent genera- 
tion in graphene at normal incidence. In bulk graphene 
sheets of the Deh or Cq v symmetry arguments the effect 
is forbidden. In order to explain its appearance in exper- 
iment, we need to take into account the graphene edges 
which locally reduce the symmetry^. 



A. Model 

The existence of edge currents induced both by linearly 
or circularly polarized radiation can be understood in the 
framework of the model depicted in Fig. [7J We will show 



that the current can be generated along the sample edge 
in a narrow channel of the width comparable to the mean 
free path I. Then, if the opposite sample edges, in the 
following the left and right edges, are non-equivalent, the 
edge currents do not compensate each other resulting in a 
non-zero net electric current. Otherwise, currents flowing 
along the opposite edges are oppositely directed and have 
the same absolute values. In the simplest model of the 
edge inequivalence, the electrons are specularly reflected 
from the right edge of the sample and diffusively scat- 
tered by the left one, as illustrated in Fig. [7J In a more 
general model, we can separate the edge reflectivity into 
two parts, specular R sp and diffuse i?dif, satisfying the 
identity R sp + R^ii = 1. In this model the net current is 
proportional to the difference Rda.L — Rdif.R of the left- 
and right-edge diffusivenesses. For distinctness, we con- 
sider the current generation near the sample left edge and 
refer to Fig. [Jj An important point to bear in mind is 
that an electron traveling in the Z-thick stripe adjacent to 
the edge randomly changes its direction of motion along 
the y axis not only in the bulk scattering processes de- 
scribed by the relaxation time t\ but also in the diffusive 
reflection at the edge. Thus, its total scattering rate is a 
sum of two terms, 



— = - + R dii Mo(-v x ) 

Ttot r d 



(28) 



where d is the distance from the electron to the left edge, 
see Fig. [JJ and O(v) is the Heaviside function. The 
factor 0(— v x ) ensures that the additional contribution 
to r t ~^ exists only for the carrier moving towards the 
edge. The dependence of the scattering time Ttot on the 
velocity direction is plotted in Fig. [Jjb). 

The x-component of electric field results in the oscil- 
lating correction to the electron velocity x-component: 



Sv x , 



where 



5vi 



erE x v 
— ( 

(1 — iujt)p 



c.c, 



(29) 



(30) 



where p is the electron momentum corresponding to its 
mean kinetic energy. The modulation of the velocity re- 
sults in the modulation of the scattering time r in agree- 
ment with Eq. (|28p ■ since electron reaches the bound- 
ary faster of slower depending on the value of eE x , see 
Fig.UJb). 

The dc current along the edge is generated by y com- 
ponent of the field and, according to Drudc theory, is 
proportional to the time-averaged product of E y Ttot' 



e z Nl-Re{E*5T\. 
p y 



(31) 



Here N is the electron concentration, Nl is the electron 
density per unit length in the vicinity of the edge, overline 
means the averaging over time, 

5v x 2 



St = Ttot - t ~ Q(v x )- 



(32) 
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is the oscillating correction to the scattering time, and, 
for simplicity, we assumed that i?diff,L ~ 1- The value 
and direction of the dc current is controlled by the rel- 
ative phase of the E y and r to t, or accordingly v x oscil- 
lations. At lot <C 1 the oscillations of E x and v x are 
in phase (no retardation) and the dc current appears at 
linear polarization of incident radiation where E x and 
E y components oscillate in phase. With an increase of 
lot the retardation effect comes into play and, in ad- 
dition to the linear photocurrent, a helicity dependent 
response appears: the oscillations of r tot appear to be 
phase-shifted with respect to the oscillations of E x and 
for the circularly polarized wave where E x and E y are 
7r/2 phase-shifted the dc current is generated. Taking 
the average one can find that 



LPGE: J" 



CPGE: J™ c 



-(E x E* + E y E* x )T 3 N, (33) 



—i(E x E* y -E y E* x yN, (34) 



at lot -C 1^ Superposition of these two current contribu- 
tion yields polarization dependence of the photocurrent 
observed in experiment. Indeed, in the experimental ge- 
ometry, where elliptically polarized radiation is obtained 
by rotation of the A/4 plate, polarization determined 
terms reduces to E X E* + E y E* = £?sin(4<^) + C cos(4y>) 
and i(E x E* — E y E x ) = Asm2ip being in agreement 
with experimental findings (see Fig. [T] and [2]). Also 
in the set-up, where linear polarization is rotated by 
the angle a, we obtain from Eq. ([33]) E X E* + E y E* = 
2B sin(2a) + 2C cos(2a) which also agrees with experi- 
ment (see Fig. [2]). 

As we mentioned above, different edges of the sample 
make different contributions to the photocurrcnts. Con- 
sider the square-shaped sample and assume that all the 
sample edges scatter electrons diffusively. Obviously, the 
currents flowing along the opposite edges of the sample 
should have opposite directions. The current topology 
depends on the radiation polarization state. In accor- 
dance with our model, the helicity driven current should 
wind in the same direction, clock-wise or counter-clock- 
wise depending on the radiation helicity, for all the sam- 
ple edges since the helicity, i(E x E* — E y E*) is preserved 
at the rotation by 90°. Hence, the distribution of helic- 
ity driven current should form a vertex, whose winding 
direction changes with the change of light helicity. How- 
ever, it is not the case for the photocurrent caused by 
linearly polarized light and described by the combination 
E X E* + E y E*. The latter changes its sign at the 90° ro- 
tation, hence, the linear photocurrcnts flow towards (or 
outwards) the same corner for adjacent edges. 



B. Microscopic theory 

The microscopic theory of this effect is developed in 
the framework of kinetic equations for the classical fre- 



quency range. The equation for the distribution function 
f(k,x,t) of electrons in the conduction band in a semi- 
infinite sample (x > 0) has the form 



dt x dx 



eE(t) Of 
h dk 



Q{f} 



(35) 



where E(t) = E e~ iult + E^e +iut is the electric field of 
the radiation, which is assumed uniform in the graphene 
plane in the geometry of normal incidence. The distri- 
bution function can be expanded in series in the electric 
field as follows 

f(k, x, t) = f (s k ) + [.h(k, x)e~ iut + c.c] + / 2 (fc, x) + ... , 

(36) 

where /o(efc) is the equilibrium distribution function, 
/i oc \E\, and f 2 <x l-^l 2 - The directed electric current 
along the structure edge is then given by 



/>oo 

J y — 4e / dx S2 f%{k, x)v y 
Jo u 



(37) 



Here, the factor 4 accounts for the spin and valley degen- 
eracy. 

As a model, we consider the simplest form of collision 
integral, 



Q{f(k,x,t)} = - 



f(k, x,t) - f (e k ) 



(38) 



and assume the diffusive elastic scattering of carriers at 
the edge, which implies that f(k,0,t) is independent of 
the direction of k for v x > and / /(fe, 0, t)v x dtpk = 0. 

To first order in the electric field, solution of Eq. ([35]) 
with the above boundary conditions has the form 



erf' 

f 1 (k,x) = --^[E -v 

1 — 1LOT 



- ( E ■ v + -E , x v ) exp 



x Q(v x ) 



(39) 



where /q = dfo(e)/de and Q(v x ) is the step function equal 
to 1 and for v x > and v x < 0, respectively. 

The equation for the second-order correction f 2 (k,x) 
to the distribution function, which gives rise to dc electric 
current, assumes the form 



df 2 (k,x) 

v x 

ox 

which yields 



—Re 

h 



E, 



dfi(k,x) 



dk 



f2(k,x) 



(40) 



/ 2 (fe, x)dx = v x T[f 2 (k, 0) - / 2 (fe, 00)] 



2er 



Rc 



dk 



dx . (41) 



By using Eqs. (1371) and (|4"Tj) we derive for the edge electric 
current 



Ju 



- 8i fl>{r 



v x v y E* d[(E ■ u)/o] 



dk 



(42) 
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d r 



(1 



) 2 dk 



{E -v + -E , x v)fo v x \ 9{v x ) 



where the above two contributions to the current stem 
from the first and second terms on the right-hand side of 
Eq. (HJ), respectively. 

In graphene, the electron kinetic energy and velocity 
are given by ek = hv\k\ and v — vk/k, respectively. 
Therefore, we obtain for the edge photocurrent contribu- 
tion from the electrons 



eW/o(0) 
y 27rft 2 [l + {lot) 2 ] 



1,7 1- {corf 
6 1+ {lot) 2 
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x {Eq, x Eq iV + E ^ v Eq x ) + — - 



+ {lot} 



(43) 



Similarly to the case of the photon drag effect the va- 
lence band holes contribution to the photocurrent can be 
obtained from Eq. (|4"3")l by the replacement e — > — e and 

/o(0) -> / (h ° le) (0) = [1 - / (0)]. Taking into account 
both electron and hole currents we obtain that the total 
electric current emerging at the sample edge is given by 
Eq. (03J), where / (0) should be replaced by [2/ (0) - 1]. 

The dependences of the edge photocurrent (|4"3"| on the 
angles ip and a are given by 



J y {<p) oc — sin Alp + C c sin 2cp, 



J y {a) oc Ci sin 2a. 



(44) 



(45) 



VI. SUMMARY 



To summarize, our observations clearly demonstrate 
that the irradiation of monolayer graphene flakes results 
in directed electric currents of different origins. In all 
our measurements a substantial contribution to the pho- 
tocurrent is driven by the photon helicity. It can be sep- 
arated into the contribution resulting from normal in- 
cidence and the one from oblique incidence. While the 
contribution from oblique incidence is related to the bulk 
material and results from the transfer of the photon an- 
gular and linear momentum to free carriers, the effect 
at normal incidence is caused by the sample edges and 
vanishes in the bulk material. Our theory describes the 
helicity driven as well as the linear-polarization driven 
photocurrents in the classical limit where the radiation 
frequency is smaller compared to the characteristic en- 
ergy of the carriers. The treatment of the general case, 
where interband transitions should be taken into account, 
is a future task. 

Note added. At the day of the submission of this 
manuscript, Entin et al. have published a pre-print^ 
devoted to the theory of the linear photon drag effect 
in graphene. Unlike in Ref. HH, where the direct inter- 
band absorption of linearly polarized light is considered, 
here we present results of mutual experimental and the- 
oretical studies of the photon drag effect in graphene un- 
der indirect intrasubband optical transitions for both lin- 
ear and circular polarizations, with emphasis on helicity- 
dependent photocurrents. 



The microscopic theory for the edge photocurrents is in 
a good agreement with experimental results. At nor- 
mal incidence the observed photoresponse is, according to 
Eq. (|4"4")) . a sum of the contributions proportional to func- 
tion sin 2ip and B sin {Aip) + C cos (iip) oc sin (4y> + £) , see 
Fig.[SJa). The additional phase £ in experimental depen- 
dences results from the arbitrary orientation of the sam- 
ple edges. Moreover, the dependence of the experimental 
data of the azimuthal angle a, Fig. [5{b), demonstrates 
sin 2a and cos 2a behavior in agreement with Eq. (|45[) . 
Hence, helicity driven and linear polarization dependent 
photocurrents in graphene arc well described at the nor- 
mal incidence by our model of edge effects. 
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